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We report on p — A, p — A, p — A and p — A correlation functions constructed in central Au-Au 
collisions at ^/s NN = 200 GeV by the STAR experiment at RHIC. The proton and lambda source 
size is inferred from the p — A and p — A correlation functions. It is found to be smaller than the pion 
source size also measured by the STAR experiment at smaller transverse masses, in agreement with 
a scenario of a strong universal collective flow. The p — A and p — A correlation functions, which are 
measured for the first time, exhibit a large anti-correlation. Annihilation channels and/or a negative 
real part of the spin-averaged scattering length must be included in the final-state interactions 
calculation to reproduce the measured correlation function. 
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I. INTRODUCTION 

Correlations amongst non-identical particles are sen- 
sitive to the space-time extent of their emitting source 
(see e.g. [H). Originally uncorrected particles pro- 
duced in nearby phase space points in the prompt emis- 
sion final state can interact through the nuclear and/or 
the Coulomb force and become correlated at time scales 
much longer than the production time. When the final- 
state interaction (FSI) is relatively well understood the 
emitting source size can be inferred from correlations 
at small relative three-velocity of the particles in their 
center-of-mass system. In relativistic heavy-ion collisions 
large particle densities are produced and the collision fire- 
ball may undergo a collective expansion (i.e. flow) 
This flow can induce space-momentum correlation so that 
particles with similar velocities come from the near-by 
regions of the source. With a strong flow at RHIC.as 
suggested by several measurements (see e.g. 0,11,13, Hi), 
the observed source sizes should be reduced relative to a 
source without flow [2!| and vary with the mass of the 
emitted particle: the heavier the particle, the smaller is 
the reduction of the collective flow effect due to the ther- 
mal motion and the smaller are the apparent source sizes 

@ . This flow effect can also be studied with p—p corre- 
lations and compared with n — n or K — K correlations. 
As compared with the p — p system, the p — A system 
gains in statistics in the region of small relative velocities 
due to the absence of repulsive Coulomb interaction Q- 

In this paper we test the hypothesis that a strong flow 
is established in Au+Au collisions at ^/s NN = 200 GeV 
by comparing the source sizes of protons and lambdas to 
that of pions. The first measurements of p — A, p — A and 
p — A correlation functions are presented. The p — A and 
p— A interaction potentials are relatively well understood, 
so we are able to infer source sizes jf| 0, H, 0, E3, El ■ The 
p — A and p — A FSI, on the other hand, are unknown. 
As such, the scattering lengths and source sizes are ex- 
tracted by fitting the data from the STAR experiment 
with the Lednicky & Lyuboshitz analytical model 
Besides constraining baryon - antibaryon potentials, this 
information determines unknown p— A andp— A annihila- 
tion cross sections which are useful to constrain heavy-ion 
cascade models fllij . 

II. DATA RECONSTRUCTION 
A. Events selection 

The analysis was carried out using the STAR detec- 
tor at RHIC (l4j . Two million Au+Au collisions have 
been analyzed with ^/s NN = 200 GeV. Because of statis- 
tics issues, only the 10% most central collisions were se- 



lected with the Zero-Degree Calorimeters and the Cen- 
tral Trigger Barrel of the STAR detector. This event 
selection procedure is explained in detailed in [f|. The 
other centrality selections didn't gave any statistically 
meaningful results. Tracking of charged particles was ac- 
complished using the STAR Time Projection Chamber 
(TPC) which covers the kinematic range of transverse 
momentum p t > 150 MeV/c, pseudo-rapidity | -^7 1 < 1.5 
and azimuthal angle < 4> < 27r. Events analyzed in 
this paper have collision vertices within ±25 cm longitu- 
dinally of the TPC center. 



B. Protons and antiprotons selection 

Protons and antiprotons are identified using their spe- 
cific energy loss (dE/dx) in the TPC gas. This selection 
limits the acceptance of particles to the transverse mo- 
mentum range of 0.4 — 1.1 GeV/c in the rapidity interval 
\y\ < 0.5. Tracks pointing to within 3 cm of the primary 
vertex are included in the primary track sample. 



C. Lambdas and antilambdas selection 

Lambdas (antilambdas) are reconstructed through the 
decay channel A — * tt~ + p (A — > n + + p) [l5j], with 
branching ratio of 64%. Pions and protons (i.e. lambda 
daughters) are selected using their specific energy loss. 
The invariant mass (Figure [1]) range of the lambda can- 
didates is 1115 ± 6 MeV/c 2 , the ±6 MeV/c 2 has been 
fixed to optimize the signal over noise ratio. The signal 
over noise ratio is equal to 86 ± 6% for lambda ((pt) = 
1.05 GeV/c) and for antilambda ((p t ) = 1.09 GeV/c) in 
the ±6 MeV/c 2 mass window and 0.3 < p t < 2.0 GeV/c. 
The correlation effect was the same within the errors with 
one sigma cut (±3 MeV/c 2 ) on the invariant mass. In 
addition the following geometrical cuts are applied. The 
distance of the closest approach (DCA) of lambda daugh- 
ters is required to be less than 0.7 cm. The DCA of the 
decay pions with respect to the primary vertex is required 
to be greater than 2.0 cm. The DCA of the reconstructed 
neutral particles to the primary vertex is required to be 
less than 0.6 cm. To avoid Kg being misidentified as 
lambdas, lambda candidates are rejected if their invari- 
ant mass is within the window 497.7^21 3 MeV/c 2 when 
the pion mass is assumed for the two daughters. Due to 
the detector acceptance and the selection criteria, the pt 
range of the lambda sample is 0.3 < pt < 2.0 GeV/c and 
|y|<1.5. 
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FIG. 1: (Color online) Invariant mass of the selected A (a) and 
A (b) background not substracted with 0.3 < pt < 2.0 GeV/c. 
The y-axis represents the number of candidates used in this 
analysis. 



D. Pairs selection 

When studying two-article FSI, the relevant variable 
is the momentum of one of the particles in the pair rest 
frame called here k* (k* = \k*\). The correlation function 
has been extracted by constructing the ratio of two dis- 
tributions. The numerator is the k* distribution of pairs 
from the same event. The denominator is the k* distribu- 
tion of pairs composed of particles from different events 
with primary vertices separated from each other by less 
than 10 cm. The ratio is formed by dividing the numera- 
tor by the denominator. Then the ratio is normalized to 
1 at high k* (k* > 0.35 GeV/c). The event mixing proce- 
dure is the same as the one used in [j| . When reconstruct- 
ing a primary lambda (antilambda) the decay proton (an- 
tiproton) points directly back to the primary vertex and 
may share some hits with a primary proton (antiproton) 
in the TPC. This phenomenon is called track merging 
and can occur while building pairs for the correlation 
function. In case of track merging, instead of counting 
two tracks with small k* only one track will be found. So 
one pair will be missed at small k* . If a lot of pairs are 
missed, the correlation function will show a hole at small 
k* because they are not found. A missed pair leads to a 
fake correlation because of the event mixing procedure. 
A pair can be missed in a real event (in the numerator). 
Such a pair may be reconstructed taking two different 
events in order to build the background (in the denom- 
inator). Thus, track merging leads to fake correlations. 
For two tracks of different momenta, or different polar 
angles, the number of shared hits may vary as a function 
of where they cross in the TPC. It could be as low as 5 



hits on the edge of the TPC. The tracker can not find 
these hits, it is linked with the finding seed. This affects 
high p t tracks because they will have more hits merged. 
To avoid such fake correlations, track merging has been 
studied for all other possible track - daughter track com- 
binations. The study of track merging for lambda daugh- 
ters and proton/antiproton tracks leads to two different 
selections criteria. The first selection requires tracks to 
share fewer than 10% of their TPC space points. The sec- 
ond selection deals with the average separation between 
primary tracks and lambda/antilambda daughter tracks. 
The track separation is calculated as an arithmetic mean 
distance between the TPC hits of the two tracks for a 
given radius. If a track crosses the whole TPC, it will be 
reconstructed with a maximum of 45 hits. Because one 
of the lambda/antilambda daughter tracks is a secondary 
track, all 45 hits of the TPC may not be available, so the 
mean is calculated from a maximum of 11 distances. In 
this paper " secondary particles" means particle from de- 
cay. As a consequence the average separation between 
protons/antiprotons and lambda/antilambda daughters 
are required to be greater than 11 cm for p — pj-, p — Pa, 
p — 7Ta, p — 7Tj, 10 cm for p — p\, 12 cm for p — p^, 17 cm 
for p — 7Tyv and p — ir-^. The first selection prevents inter- 
ference between opposite sign tracks that, even though 
their average separation is large, can cross each other in 
the TPC. When the trajectories of the particles cross, 
space points can be assigned to the wrong track during 
reconstruction. In some cases this can lead to a failure 
to reconstruct one of the tracks. For this reason the val- 
ues for the minimum average separation are larger for 
opposite sign pairs. 



III. PURITY 

A. Definition 

Impurities in the sample of protons and lambdas will 
reduce the observed p — A correlation strength. In the 
case of lambdas, fake lambda candidates (from combina- 
torial background) and secondary lambdas (e.g. from S° 
decays) are the two main sources of impurity. The sam- 
ple of protons is contaminated by other charged tracks 
falsely identified as protons and by protons from weak 
decays (feed-down). In order to correct the observed cor- 
relations for misidentification and feed-down we estimate 
the particle purity for p, p, A and A as a function of trans- 
verse momentum (pt)' 

Particle Purity(p t ) = Pid(p t ) x Fp(p t ), (1) 

where Pid is the probability a candidate was correctly 
identified and Fp is the fraction of the candidates that 
were primary particles. The final correction depends on 
the product of the particle purity for both particles (i.e. 
the Pair Purity). The pair inpurity is corrected for in 
constructing the correlation function in k* . 
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The feed-down estimations have been done for p, p, A 
and A and are summarised in Table |U Combined results 
from STAR [ 3 ,1^ . JlR Il7j and predictions from a thermal 
model fl8l . fl9L |20|] have been used. The approximations 
introduced by estimating the purity are the major source 
of systematic uncertainties on the extracted values of FSI 
parameters and source radii discussed below. 



B. Protons and antiprotons purities 



Particle 


Identification 


Fraction Primary 


P 


76 ± 7% 


52 ± 4% 


P 


74 ± 7% 


48 ± 4% 


A 


86 ± 6% 


45 ± 4% 


A 


86 ± 6% 


45 ± 4% 



TABLE I: Summary of the particle purity due to identifi- 
cations and weak-decay contamination. Values are averaged 
over the transverse momentum without taking into account 
the transverse momentum dependence for k* < 0.2 GeV/c. 



The identification probabilities have been estimated 
for charged particles; they are also given in Table HI 
A track can be identified as a pion, a kaon, a pro- 
ton/ antiproton and an electron/positron with a certain 
probability using the information about the energy loss 
dE/dx [21(. Identified protons (antiprotons) from the se- 
lected sample have a mean identification probability of 
76 ± 7% (74 ± 7%). This mean identification probability 
and the corresponding p t are calculated over all selected 
tracks considered as protons (antiprotons). The calcu- 
lated feed-down leads to a mean estimated fraction of 
primary protons of 52% (with a mean transverse velocity 
{Pt) = (pt/j)/m) = 0.58, m is the mass of the parti- 
cle). Most of the secondary protons come from lambda 
decays (primary A and A from £°, S° and S - ) and con- 
stitute 36% of the protons used to construct the correla- 
tion function. Other sources of contamination for protons 
are products of S + decays and interactions of pions with 
detector materials which represent respectively 10% and 
2% of the sample. The feed-down study for antiprotons 
{{Pt) — 0.60), leads to an estimated fraction of primary 
antiprotons of 48%. Most of the secondary antiprotons 
come from antilambda decay (primary A and A from S°, 
S° and S~) and constitute 39% of the antiprotons used 
to construct the correlation function. Antiprotons from 
£+ decays are another major source of contamination; 
they make up 13% of the antiproton sample. 



C. Lambdas and antilambdas purities 

For lambdas and antilambdas the probability of 
misidentification corresponds to background estimation 
under the mass peak. The corresponding identification 
probabilities are practically independent of pt and equal 
to 86 ± 6% for both lambdas and antilambdas, respec- 
tively. The sample of lambdas (antilambdas) includes 
s eco ndary particles such as decay products of H~, S°, S° 
(H _ ,S°,S°). The fractions of primary lambdas {{Pt) = 
0.68) and primary antilambdas {{Pt) = 0.70) have been 
estimated at 45%. 



D. Pairs purities 

The pair purity plays a crucial role in the correlation 
study. The estimated value of the mean pair purity for 



p — A, p — A, p — A and p — A systems is A = 17.5 ± 
2.5% after taking into account the transverse momentum 
dependence. Without taking into account the transverse 
momentum dependence the estimated purities differ by 
2% (Tab. [fl. 



IV. CORRECTIONS 
A. Purity 

Because the contamination reduces the correlation 
strength, raw data have been corrected for purity using 
the relation: 



c 



corr 

measured 



-^measure 



d{k*) 



1 



PairPurity 



1, 



(2) 



where PairPurity is the product of the purities for the 
two particles and C^ sured {k*) and C m 

easured 

(fc*) are 

respectively the corrected and measured correlation func- 
tions. Eq. ^ assumes that misidentified and weak decay 
protons (antiprotons) are uncorrelated with lambdas and 
antilambdas. This assumption is justified for misidenti- 
fied protons (antiprotons) since the eventual pion or kaon 
correlation at small k* is washed out after the wrong mass 
assignment. Combinatoric background reconstructed as 
A and A also leads to uncorrelated pairs. On the other 
hand weak decay products may keep a residual correla- 
tion from their parents. This assumption will be revisited 
when extracting source sizes and scattering lengths from 
the correlation functions. 



B. Momentum resolution 

The effects of momentum resolution have been stud- 
ied using mixed pairs by calculating the weig ht with the 
Lednicky & Lyuboshitz analytical model [12| ■ It appears 
that compared with statistical and systematic errors, the 
impact of the momentum resolution effect is negligible. 
Indeed, the momentum resolution effect leads to about 
1% variation of the apparent source radius. Nevertheless, 
correlation functions have been corrected for the momen- 
tum resolution using the following formula: 



Ctrue (^ 



/~icorr 
measure 



x CTh-not — smeared 



Crh- 



smearei 



>{k* 



, (3) 
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where Ct rue (fc*) represents the corrected correlation 

function, CTh-not-smeared(k*)/C>Th-smeared(k*) IS the 

correction factor; CTh-not-smeared(k*) is calculated 
without taking into account the effect of momentum reso- 
lution and CTh-smeared{k*) includes this effect. The shift 
due to the momentum resolution is studied using simu- 
lated tracks introduced into real events. This shift is ap- 
plied to momenta in order to calculate Crh-smearedik*). 
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FIG. 2: (Color online) The purity and momentum-resolution 
corrected correlation functions Ctme(k*) for p — A, p — A (a), 
p — A, p — A (b). Curves correspond to fits done using the 
Lednicky & Lyuboshitz analytical model [l2l ]. 



RESULTS 



Correlation functions 



In FigTJ] (a) the corrected p — A and p — A correla- 
tion functions are shown. They are close to each other, 
within error bars, showing a pair excess at small k* 
(0 < k* < 0.1 GeV/c) indicating an attractive potential 
between (anti)proton and (anti)lambda. FigJ5] (b) shows 
the corrected p — A and p — A correlation functions mea- 
sured for the first time. They are below unity in a wide 
k* range < k* < 0.25 GeV/c consistent with positive 
imaginary parts of the s-wave scattering lengths (due to 
the open annihilation channels) and a negative real part 
of the spin averaged s-wave scattering length. In Figure 
[3] the combined (p - A) © (p - A) and (p - A) © (p - A) 
correlation functions are presented. The symbol means 
that numerators and denominators of the systems have 
been added to build the combined correlation functions. 
In both Figures, curves correspond to a fit carried out 
with the Lednicky & Lyuboshitz analytical model [l2^ . 



B. Lednicky & Lyuboshitz analytical model 

This model relates the two-particle correlation func- 
tions with source sizes and scattering amplitudes [H, H2| • 
As usual, similar to the Fermi factor in the theory of [3- 
decay, the correlation function (C(k*)) is calculated as 
the square of the wave function (^ s ) averaged over the 
total spin S and over the distribution of relative distance 
(r*) of particle emission points in the pair rest frame: 



(4) 



It should be noted that the two particles are generally 
produced at non-equal times in their center-of-mass sys- 
tem and that the wave function in Eq. (01 should be 
substituted by the Bethe-Salpeter amplitude. The latter 
depends on both space (r*) and time (t*) separation of 
the emission points in the pair rest frame and at small \t* \ 
coincides with the wave function ^> s up to a correction 
0(\t* /mr* 2 \), where m is the mass of the lighter particle. 
It can be shown that Eq. ((4]) is usually valid better than 
to few percent even for particles as light as pions [ll, HH . 
The wave function VP represents a stationary solution of 
the scattering problem having at large distances r* the 
asymptotic form of a superposition of the plane and out- 
going spherical waves. It is approximated by the solution 
outside the range of the strong interaction potential tak- 
ing into account, at the considered small k* values, the 
s-wave part of the scattered wave only: 



(5) 



with the effective range approximation for the s-wave 
scattering amplitude: 



f(k* 



Jo L 



(6) 



where f$ is the scattering length and d,Q is the effective 
radius for a given total spin S=l or S=0, i.e. for a triplet 
(t) or singlet (s) state respectively. We assume that par- 
ticles are produced unpolarized, i.e. po = 1/4 of the pairs 
are in the singlet state and p\ = 3/4 are in the triplet 
state. Then, assuming a Gaussian distribution in r* , 



d 3 N/d 3 r* 



(7) 



where ro can be considered as the effective radius of the 
source, the correlation function can be calculated analyt- 
ically [13: 



c{k*) = i + J2ps 

s 

2^f s (k*) 



f S (k*) 



ro 



1 - 



2y/jrr 



Fi(Qro) - J y > F 2 (Qr ) 



(8) 



Vnr r 

where Fi(z) = f*dxe x2 ~ z2 /z and F 2 (z) — (1 — 

e _z )/z. The leading correction to the correlation func- 
tion 0{\f^\ 2 dQ /tq) is introduced in Eq. (jHJl to account 
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for the deviation of the solution ([5]) from the true wave 
function inside the range of the strong interaction poten- 
tial. 
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FIG. 3: (Color online) (p - A) © (p - A) and (p - A) © (p - A) 
combined correlation functions. Correlation functions are cor- 
rected for purity and momentum resolution. Curves corre- 
spond to fits done using the Lednicky & Lyuboshitz analytical 
model El- 



C. FSI parameters and source sizes 

The p — A and p — A interactionpotentials are rela- 
tively well understood [1, 0, H, H, ES E3 , which allows us 
to extract the source radius tq from the fit. The best fits 
are compared with the separate p — A and with the p — A 
correlation functions in Figure [5] (a), and the combined 
one in Figure [31 The scattering lengths (/q — 2.88 fm, 
fl = 1.66 fm) and effective radii (dg = 2.92 fm, 4_= 3.78 
fm) from 0] have been used for the p — A, p — A corre- 
lation functions. The systematic errors on the radius 
introduced by the uncertainties on the scattering lengths 
have been estimated to be 0.2 fm assuming spin averaged 
FSI parameters with 5% uncertainty. The fit results are 
summarized in Table [II] The three errors are, from left 
to right, the statistical errors and the systematic errors 
introduced by the uncertainty on the purity correction 
and on the scattering length for p — A and p — A systems 
and on the uncertainty in the model for p — A and p — A 
systems.. One parameter is free while fitting the p — A, 
the p — A and the combined correlation functions. Three 
parameters are free while fitting the p — A, the p — A and 
the combined correlation functions. Statistical errors on 
the radii are larger for the p — A, the p — A and for the 
combined correlation functions than for the correspond- 
ing baryon-antibaryon ones. The p — A, the p — A and the 
combined correlation functions have a large width and in- 



volve more statistics in the fit of the correlated fc*-region 
as compared with the p — A, the p — A and the combined 
correlation functions. 

The extracted source radii are close to the values (3 — 
4 fm) obtained in measurements performed by the NA49 
(SPS) collaboration in Pb+Pb collisions at 158 AGeV @ 
and by the E895 (AGS) experiment in Au+Au collisions 
at 4, 6, and 8 AGeV @, [Hj]. This confirms that the 
particle emitting source size does not change significantly 
with beam energy; a result also obtained by studying 
two-pion correlations. 

The p — A and p — A scattering lengths have never 
been measured before. Hence, they have to be included 
as free parameters in the fit to the experimental corre- 
lation functions. In order to limit the number of free 
parameters, the following assumptions are made: i) the 
spin dependence is neglected, f s = /* = / and ii) the 
effective radius (do) is set to zero. An extra parameter 
Im/o > is added, taking into account the annihilation 
channels. 

The best fits are compared with the separate p — A and 
p — A correlation functions in Figure [2] (b) , and with the 
combined one in Figure O The fitted spin- averaged scat- 
tering lengths for the combined p — A and p — A systems 
are compared with measurements for the p — p system 
[H, HI, III, H3] in Figure H The imaginary part of the fit- 
ted scattering length is in agreement with the p— p results 
while the real part is more negative. The error contour 
represents the statistical errors. The systematic error 
due to the uncertainty on pair purity is investigated by 
comparing the best estimated fc*-dependent purity with 
k* -independent purity corrections A. The fc*-dependent 
purity correction tends to decrease the size of the er- 
ror contour (the curve labeled "Corrected" ) as compared 
with the one for the constant purity (the curve labeled 
"A = 17.5%") and shifts both the real and imaginary 
parts of the scattering length in the direction of the p — p 
values. 

The radii extracted from the fits to the separate and 
the combined p — A and p — A correlation functions are 
summarized in Table [TTJ The errors include from left to 
right, statistical errors, systematic errors due to purity, 
and systematic errors estimated from varying model pa- 
rameters. The error on the radius parameter due to the 
uncertainties on the model is estimated to be 0.3 fm. This 
error is estimated by fixing the spin averaged scattering 
lengths and by extracting the radius and the effective ra- 
dius. For the moment, we do not have any tool to extract 
the radius uncertainty related to the neglected p-wave 
contribution. A larger radius implies a correlation over 
a smaller fc*-range than seen in the data, which cannot 
be recovered by increasing the magnitude of scattering 
lengths. However the radii extracted from the p — A (and 
p— A) and thep— A (andp— A) are significantly different. 
The error bars accounting for all statistic and systematic 
contributions barely overlap. 
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FIG. 4: (Color online) The combined (p — A) © (p — A) spin- 
averaged s-wave scattering length compared with the previous 
measurements for the p — p system [24|, HE EE E3- The 
curves show the one standard deviation contours. Note that 
for (p — A) © (p — A) only, one should read 0.1973 x Im/o 
instead of Im/o on the x-axis. 



Exp. 


System 


r (fm) 


STAR 


p- A 


2.97 ± OMto.H ± 0.2 


STAR 


p - A 


3.24 ±0.59lo;i4 ±0.2 


STAR 


p- Affip- A 


3.09 ± 0.301^25 ± 0.2 


STAR 


p - A 


1.56 ± 0.08±";"i ± 0.3 


STAR 


p- A 


1.41 ±0.10 ±0.11 ±0.3 


STAR 


p — A ©p - A 


1.50 ± 0.05±g;„ ± 0.3 


NA49 


p — A (A = 0.33 fixed) 


3.8 ±0.33 


NA49 


p-A (A = 0.17 ± 0.11 free) 


2.9 ±0.7 


E895 


p- A (A = 0.5 ±0.2 free) 


4.5 ±0.7 



TABLE II: Comparison of the radius of the source of particles 
for p — A, p— A, p— A, p— A and combined systems. For STAR, 
the three errors are, from left to right, the statistical errors 
and the systematic errors introduced by the uncertainty on 
the purity correction and on the scattering length for p — A 
and p — A systems and on the uncertainty in the model for 
p — A and p — A systems. For NA49 @] and E895 0, the A 
parameter represents the pair purity. 



VI. DISCUSSION 

The difference in radii between p — A (p — A) and p — A 
(p — A) is unexpected. Indeed, it would imply a novel 
dynamical space-momentum correlation between proton 
(p) and A (A). Strong space- momentum correlations are 
exhibited in Au±Au collisions at RHIC. These are un- 
derstood to arise from the collective flow of massive par- 
ticles [E H^|. This effect, however, would not lead to a 



difference between the source size measured from p — A 
and p — A correlations. In Figure 4, the source sizes 
from proton-Lambda correlations and pion-pion correla- 
tions are plotted as a function of the mean of the par- 
ticles' transverse masses. The decrease of the source 
size with increasing mean transverse mass is in quali- 
tative agreement with expectations from collective flow 
[29l | . The curve in FigfS] represents an arbitrarily normal- 
ized (mt)" 1 / 2 dependence. This dependence is expected 
within some hydrodynamics- motivated models [301 ] - The 
data are in reasonable agreement with this expectation. 
In addition, a possible difference between radii pointed 
by data may imply that baryon - antibaryon pairs are 
produced close in space, a dynamic correlation that is 
not in baryon - baryon pairs. 
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FIG. 5: (Color online) Pion source size [f| compared with 
proton and lambda source sizes as a function of the mean 
transverse mass ({m t )). The curve shows the (m t ) -1 '' 2 de- 
pendence with an arbitrary normalization. 

While a novel space-momentum correlation between 
proton and A cannot be ruled out, the difference between 
the radii extracted from p — A and p — A correlation func- 
tions may come from an imperfect treatment of the purity 
correction. Indeed, we have assumed that any pairs that 
are not composed of two primary particles are not corre- 
lated. However, Table HTT1 shows that a number of such 
jairs may carry a residual correlation from their parents 
3~lj ] - For example, the interaction between a primary pro- 
ton and a E° may not be completely washed out when 
constructing the p — A correlation function with the A 
being the S° daughter as the A carries most of the mo- 
mentum of its parent. This effect was found to be on 
the order of 10% [Hj]. However, none of the interactions 
between the pairs listed in Table IIIII are known. We are 
thus unable to perform any reliable correction or error 
estimate. At that stage, we show the p — A and p — A 
correlation function corrected with the best estimate of 



9 



the purity assuming no residual correlations. We extract 
source radii and scattering length parameters acknowl- 
edging that the values may be biased by the presence of 
residual correlations. 



Pairs 


Fractions 


Pprim ^prim 


15% 


PA A-prim 


10% 


??E + A-prim 


3% 


Pprim ~~ -A-E 


11% 


p A - A s0 


7% 




2% 


Pprim -A-H 


9% 


Pa -A H 


5% 


Pe+ ~ A s 


2% 


Pprim Pprim 


7% 



TABLE III: Summary of the main fractions of pairs contain- 
ing particles from particle decays included in p — A, p — A, 
p — A, p — A correlation functions assuming the absence of 
residual correlations. Ah, are A (A) decay products of H~, S° 
(H~,H°), A E o, are A (A) decay products of S° (E°), pa are 
p (p) decay products of A (A), p E + are p (p) decay products 
of £+ (£+), Aprim and Pprim represent primary A (A) and 
p (p). The remaining 29% represents misidentified p (p) and 
reconstructed fake A (A). 



VII. CONCLUSION 

Constructing p — A, p — A, p — A, p — A, we have gath- 
ered information about the space-time features of baryon 



and antibaryon emission and about the interaction in 
p — A and p — A systems. The source radii extracted 
from p — A and p — A corrrelation function agree with 
the flow expectation. The radii extracted from p — A 
and p — A are significantly smaller. Final-state interac- 
tions parameters, such as spin averaged s-wave scattering 
length, have been extracted from p — A and p — A cor- 
relation functions. The real part of the scattering length 
appears to be negative while the imaginary part is posi- 
tive, the latter being required by the unitarity due to the 
open annihilation channels. These results demonstrate 
that correlation measurements can be used to study the 
two-particle strong interaction for particle combinations 
that are difficult to access by other means, including tra- 
ditional scattering experiments. 



Acknowledgments 

We thank the RHIC Operations Group and RCF at 
BNL, and the NERSC Center at LBNL for their sup- 
port. This work was supported in part by the HENP 
Divisions of the Office of Science of the U.S. DOE; the 
U.S. NSF; the BMBF of Germany; IN2P3, RA, RPL, 
and EMN of France; EPSRC of the United Kingdom; 
FAPESP of Brazil; the Russian Ministry of Science and 
Technology; the Ministry of Education and the NNSFC 
of China; 1RP and GA of the Czech Republic, FOM of 
the Netherlands, DAE, DST, and CSIR of the Govern- 
ment of India; Swiss NSF; the Polish State Committee 
for Scientific Research; STAA of Slovakia, and the Korea 
Sci. & Eng. Foundation. 



[1] R. Lednicky, V.L. Lyuboshitz, B. Erazmus and 

D. Nouais, Phys. Lett. B373, 30 (1996). 
[2] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 112301 

(2004). 

[3] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 182301 
(2004). 

[4] J. Adams et al. (STAR), Phys. Rev. C72, 014904 (2005). 
[5] J. Adams et al. (STAR), Phys. Rev. C71, 044906 (2005). 
[6] C. Blume et al. (NA49), Nucl. Phys. A715, 55 (2003). 
[7] F. Wang and S. Pratt, Phys. Rev. Lett. 83, 3138 (1999). 
[8] P. Chung et al. (E895), Phys.Rev. Lett. 91, 162301 
(2003). 

[9] Th.A. Rijken, V.G.J. Stoks and Y. Yamamoto, Phys. 

Rev. C59, 21 (1999). 
[10] P.M.M. Maessen, Th.A. Rijken and J.J. de Swart, Phys. 

Rev. C40, 2226 (1989). 
[11] M.M. Nagels, Th.A. Rijken and J.J. de Swart, Phys. Rev. 

D20, 1633 (1979). 
[12] R. Lednicky, V.L. Lyuboshitz, Sov. J. Nucl. Phys., 35, 

770 (1982); Proc. CORINNE 90 Nantes, France (ed. Ar- 

douin, World Scientific), 42 (1990). 
[13] G.J. Wang, G. Welke, R. Bellwied and C. Pruneau, nucl- 



th/9807036. 

[14] K.H. Ackermann et al. (STAR), Nucl. Phys. A661, 681c 
(1999); K.H. Ackermann et al. (STAR), Nucl. Instrum. 
Meth. A499, 624 (2003). 

[15] J. Adams et al. (STAR), Phys. Lett. B567, 167 (2003). 

[16] C. Adler et al. (STAR), Phys. Rev. Lett. 89, 092301 
(2002). 

[17] M.A.C. Lamont, Ph.D. Thesis University of Birmingham 
(2002). 

[18] P. Braun-Munzinger and J. Stachel, Nucl. Phys. A606, 
320 (1996). 

[19] P. Braun-Munzinger, I. Heppe, and J. Stachel, Phys. 

Lett. B465, 15 (1999). 
[20] P. Braun-Munzinger, D. Magestro, K. Redlich and 

J. Stachel, Phys. Lett. B518, 41 (2001). 
[21] C. Adler et al. (STAR), Phys. Rev. Lett. 87, 182301 

(2001). 

[22] R. Lednicky, Proc. of the International Workshop on the 

Physics of the Quark-Gluon Plasma, Palaiseau, (2001). 
[23] R. Lednicky, nucl-th/0501065. 

[24] I.L. Grach, B.O. Kerbikov and Yu.A. Simonov, Sov. J. 
Nucl. Phys. 48, 609 (1988). 



10 



[25] C.J. Batty, Rep. Prog. Phys. 52, 1165 (1989). 

[26] H.J. Pirner, B.O. Kerbikov and J. Mahalanabis, Z. Phys. 

A338, 111 (1991). 
[27] E. Klcmpt, F. Bradamante, A. Martin and J.M. Richard, 

Phys. Rep. 368, 119 (2002). 
[28] J. Adams et al. (STAR), Phys. Rev. Lett. 91, 262302 

(2003). 

[29] F. Retiere and MA. Lisa, Phys. Rev. C70, 044907 



(2004). 

[30] B. Tomasik, U.A. Wiedmann, and U. Heinz, Nucl. Phys. 

A663, 753 (2000). 
[31] A.V. Stavinsky et al., Rapport interne Subatech-2003-03 

(2003). 

[32] F. Wang, Phys. Rev. C60, 067901 (1999). 



